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We theoretically study the kinetics of photoinduced reordering triggered by linearly 
polarized (LP) reorienting light in thin azo-dye films that were initially illuminated 
with LP ultraviolet (UV) pumping beam. The process of reordering is treated as a 
rotational diffusion of molecules in the light intensity-dependent mean-field poten- 
tial. The two dimensional (2D) diffusion model which is based on the free energy 
rotational Fokker-Planck equation and describes the regime of in-plane reorientation 
is generalized to analyze the dynamics of the azo-dye order parameter tensor at vary- 
ing polarization azimuth of the reorienting light. It is found that, in the photosteady 
state, the intensity of LP reorienting light determines the scalar order parameter 
(the largest eigenvalue of the order parameter tensor) , whereas the steady state ori- 
entation of the corresponding eigenvector (the in-plane principal axis) depends solely 
on the polarization azimuth. We show that, under certain conditions, reorientation 
takes place only if the reorienting light intensity exceeds its critical value. Such 
threshold behavior is predicted to occur in the bistability region provided that the 
initial principal axis lies in the polarization plane of reorienting light. The model 
is used to interpret the experimental data on the light-induced azimuthal gliding of 
liquid-crystal easy axis on photoaligned azo-dye substrates. 

PACS numbers: 61.30.Hn, 42.70.Gi 

Keywords: nematic liquid crystal; easy axis gliding; photo-alignment 



I. INTRODUCTION 

Linearly polarized (LP) ultraviolet (UV) or visible light is known to have a profound effect 
on the physical properties of some photosensitive materials such as compounds containing 
azobenzene and its derivatives [lj. In particular, such materials may become dichroic and 
birefringent under the action of LPUV light. Over the past few decades this phenomenon - 
the so-called effect of photoinduced optical anisotropy (POA) — which has a long history 
dating back almost a century to the original paper by Weigert ^ has been attracted much 
attention for both technological and fundamental reasons. POA is of primary importance 
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in the development of tools dealing with the light-controlled anisotropy and the materials 
that exhibit POA are very promising for use in many photonic applications 

POA also lies at the heart of the photoalignment (PA) technique which is employed in 
the manufacturing process of liquid crystal displays for fabricating high quality aligning 
substrates 0, @| ■ This method avoids the drawbacks of the traditional mechanical surface 
treatment and uses LPUV light to induce anisotropy of the angular distribution of molecules 
in a photosensitive film jsl, [9]. 

When the irradiated layer is brought in contact with a liquid crystal (LC), the surface 
ordering originated from the photoinduced anisotropy manifests itself in the anisotropic part 
of the surface tension known as the anchoring energy and determines the anchoring charac- 
teristics such as the anchoring strengths and the easy axis, n e , that defines the direction of 
preferential orientation of LC molecules at the surface. In this way, the light can be used as 
a means to control the anchoring properties of photosensitive materials. 

This is the light-induced anisotropy of the orientational distribution of molecules that 
reveals itself in POA and can be described as the photoinduced orientational ordering in 
photosensitive materials. Such ordering, though not being understood very well, can gener- 
ally occur by a variety of photochemically induced processes. These typically may involve 
such transformations as photoisomerization, crosslinking, photodimerization and photode- 
composition (a recent review can be found in Refs. [g, 10]). 

The mechanism of photoinduced ordering underlying POA and related photoaligning 
properties therefore cannot be universal. Rather they crucially depend on the material in 
question and on a number of additional factors such as irradiation conditions, surface inter- 
actions etc. These factors combined with the action of light may result in different regimes 
of the photoinduced ordering kinetics leading to the formation of various photoinduced ori- 
entational structures. 

The photoalignment has been studied in a number of polymer systems including dye 
doped polymer layers jll|, 12], cinnamate polymer derivatives (131-161" and side chain poly- 
mers containing chemically linked azochromophores (azopolymers) |17 - 20 1 . POA observed 
in similar polymers (H, 0, 17-19, 21, 22| was found to be long term stable so that the pho- 
toinduced anisotropy does not disappear after switching off the irradiation. 

Most of the above mentioned polymer systems represent azocompounds that exhibit 
POA driven by the trans- cis photoisomerization. In theoretical models (22l-l28l|. they are 
treated as ensembles of the stable trans isomers characterized by elongated rod-like molecular 
conformation and the bent banana-like shaped cis isomers. 

The mechanism of photoisomerization implies that the key processes behind the pho- 
toinduced orientational ordering of azo-dye molecules are photochemically induced trans- cis 
isomerization and subsequent thermal and/or photochemical cis-trans back isomerization 
of azobenzene chromophores. Owing to pronounced absorption dichroism of photoactive 
groups, the photoisomerization rate strongly depends on orientation of the azo-dye molecules 
relative to the polarization vector of the pumping light, ~Eiuv- 

When the cis isomers are short-living, the cis state becomes temporary populated during 
photoisomerization but reacts immediately back to the stable trans isomeric form. The 
trans- cis-trans isomerization cycles are accompanied by rotations of the azo-dye molecules 
that tend to minimize the absorption and become oriented along directions normal to the 
polarization vector of the exciting light ~Eijjv (the molecules with the optical transition dipole 
moment oriented perpendicular to ~Eijjv are almost inactive). Non-photoactive groups may 
then undergo reorientation due to cooperative motion [l8|, 25, 
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From the above it might be concluded that the photoinduced orientational structures 
will be uniaxially anisotropic with the optical axis directed along the polarization vector, 
Eyy. Experimentally, this is, however, not the case. For example, constraints imposed by 
a medium may suppress out-of-plane reorientation of the azobenzene chromophores giving 
rise to the structures with strongly preferred in-plane alignment (22|. Another symmetry 
breaking effect induced by polymeric environment is that the photoinduced orientational 
structures can be biaxial [2l|, |22j, I27L l3CH32| (a recent review of medium effects on photo- 
chemical processes can be found in [33]). 

It was recently found that, similar to the polymer systems, the long-term stable POA in 
the films containing photochemically stable azo-dye structures (azobenzene sulfuric dyes) is 



characterized by the biaxial photoinduced structures with favored in-plane alignment |34j . 135 



Unlike azopolymers, photochromism in these films is extremely weak so that it is very 
difficult to unambiguously detect the presence of a noticeable fraction of cis isomers. 

Over the last decade the photoaligning properties of such azo-dye (SD1) films has been 
the subject of intense studies 0, [36l-l38|. It was found that, owing to high degree of the 
photoinduced ordering, these films used as aligning substrates are characterized by the 
anchoring energy strengths comparable to the rubbed polyimide films. For these materials, 
the voltage holding ratio and thermal stability of the alignment turned out to be high. The 
azo-dye films are thus promising materials for applications in liquid crystal devices. 

The anchoring characteristics of the azo-dye films such as the polar and azimuthal an- 
choring energies are strongly influenced by the photoinduced ordering. In particular, the 
easy axis is dictated by the polarization of the pumping LPUV light, whereas the azimuthal 
and polar anchoring strengths may depend on a number of the governing parameters such 
as the wavelength and the irradiation dose j38|. 

In a LC cell with the initially irradiated layer, the easy axis may rotate under the action 
of secondary irradiation with reorienting LP light which polarization differs from the one 
used to prepare the aligning layer. This effect — the so-called light-induced reorientation of 
the easy axis — is governed by the photoinduced reordering of azo-dye molecules triggered 
by the reorienting light and may be of considerable interest for applications such as LC 



rewritable devices [39(. Note that slow reorientation of the easy axis (the light-induced easy 
axis gliding) on the photosensitive layers prepared using the PA technique was originally 
observed on poly-(vinyl)-alcohol (PVA) coatings with embedded azo-dye molecules [40|. For 
the azo-dye films, similar results were reported in Ref. (4ll |. 

In this paper the kinetics of photoinduced reordering underlying the light-induced easy 
axis reorientation on the azo-dye films will be our primary concern. In particular, we apply 
a generalized version of the theoretical analysis presented in Ref. [35[ to examine how the 
initial ordering of azo-dye molecules and the characteristics of reorienting light affect the 
regime of kinetics and the photosteady states. 

There are two phenomenological models formulated in Ref. 35| by using different theo- 
retical approaches: (a) the two-state model is based on the above discussed mechanism of 
photoisomerization with short living excited cis state which is characterized by weak pho- 
tochromism and negligibly small fraction of cis molecules that rapidly decays after switching 
off irradiation; and (b) the diffusion model describes the light-induced reorientation of azo- 
dye molecules as rotational Brownian motion governed by the light intensity-dependent 
mean-field potential. It turned out that predictions of the models are quite similar and they 
both were successfully used to interpret the experimental data. This equivalence points to 
the fact that, in the limiting case of short-living cis states, the photoisomerization cycles 
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with orientation dependent isomerization rates may generally be viewed as a random walk 
on the sphere described by the Fokker-Planck equation for rotational diffusion. For our 
purposes, the latter will be conveniently employed to perform the theoretical analysis and 
to model the kinetics of reordering. 

The layout of the paper is as follows. In Sec. [TT1 we briefly recapitulate the theory j35[ 
introducing the diffusion models where the process of photoinduced reordering is treated 
as a rotational diffusion in the effective mean-field potential. Then, in Sec. Ill B\ we dis- 
cuss the regime of purely in-plane reorientation and analyze the general properties of the 
corresponding two dimensional (2D) diffusion model. In Sec. II III these results are used to 
study bifurcations of the photosteady states and the kinetics of photoinduced reordering 
depending on the initial ordering, the light intensity, Ijjy and the polarization azimuth, a p , 
that characterizes orientation of the polarization vector of reorienting LPUV light, ~Eiuv- 
In Sec. [IV] we apply the model to interpret the experimental data on the light-induced az- 
imuthal gliding of the easy axis on the photoaligned SD1 substrate measured as part of 
the investigation into effects of the polarization azimuth in the dynamics of the electrically 
assisted light-induced gliding [421 ] . Finally, in Sec. |V]we draw the results together and make 
some concluding remarks. 





Figure 1: (Color online) Frame of reference: the z axis is normal to the substrate and the 
polarization vector of the activating light makes the angle a p with the y axis. 



II. MODEL 



In this section our starting point is the approach of Ref. [35j in which the light-induced 
reorientation of azo-dye molecules is described as rotational Brownian motion governed by 
the effective mean-field potential. 
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Following Ref. [35], we assume that azo-dye molecules are cylindrically symmetric, so that 
orientation of a molecule can be specified by the unit vector, u = (sin 9 cos 0, sin 9 sin 0, cos 9), 
directed along the long molecular axis. Quadrupolar orientational ordering of such molecules 



is characterized using the traceless symmetric second-rank tensor [43] 

Q(u) = (3u®u-I 3 )/2, (1) 

where I n is the n x n identity matrix. The dyadic (CQ) can be averaged over orientation 
of molecules with the normalized angular distribution function f(u,t) to yield the order 
parameter tensor S 

(Q) = J Q(u)/(u,t)du = S, (2) 

where du = sin 9d9d(p. 

A. Rotational mean-field Fokker-Planck equation 

The key assumption taken in the diffusion models formulated in Ref. [35| is that the 
angular distribution function, /(u, t), of rodlike azo-dye molecules must satisfy the rotational 
free energy Fokker-Planck (FP) equation 

ai/ .-^< r '{ /A ^fl},-L.D Iot {/.L^m}, (3) 

where D^™^ is the rotational diffusion tensor and L is the angular momentum operator that 
can be expressed in terms of of the azimuthal and zenithal (polar) angles, and 9, as follows 

iL = [rxV]=e^ e - [sin 9]~ 1 e fy, (4) 

where e# = (cos 9 cos 0, cos 9 sin 0, — sin 9) and = ( — sin 0, cos 0, 0). 

When the effective free energy functional is a sum of the effective internal energy, U[f], 
and the Boltzmann entropy, (In/), F[f] = U[f] + (In/), the free energy FP (j3J) can be 
recast into the form of the mean-field FP equation 

d t f = -L--D rot {Lf + f-LV}, V=%- (5) 

describing the rotational diffusion governed by the effective mean-field potential, V. 

1. Effective mean- field potential 

In the lowest order approximation based on the truncated expansion for the internal 
energy functional U[f] retaining one-particle (linear) and two-particle (quadratic) terms, 
the effective potential, V, is given by 

V(n) = ^-y = E7i(u) + J U 2 (u, u')/(u') du', (6) 
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where E/i(u) is the external field potential and the last term proportional to the symmetrized 
two-particle kernel, ^(u, u'), represents the contribution coming from the intermolecular 
interactions. 

The one-particle part of the effective potential (jBJ) can be written as a sum of the light- 
induced contribution 



U I (u)=u I E* uv -Q(u)-E 



UV-, 



(7) 



where an asterisk indicates complex conjugation, that comes from the interaction of azo- 
molecules with the reorienting UV light and the surface-induced potential 



U s (u) =u s z- Q(u) • z = u s Q zz (u) 



(8) 



that takes into account conditions at the bounding surfaces of the azo-dye layer. So, for the 
two-particle interaction taken in the Maier-Saupe form 

U 2 (ui, u 2 ) = u 2 Tr[Q(ui) • Q(u 2 )] = u 2 Qij{yn)Qij{vL 2 ), (9) 

we obtain the following expression for the effective potential of azo-dye molecules: 

V(u) = uJuv e p ■ Q(u) • e p + u s Q zz (u) + u 2 Sij Qy(u), (10) 

where luy is the intensity of the reorienting UV light which is assumed to be linearly 
polarized along the unit polarization vector e p 



R-r 



a, 



Rrot(«?; 



cos a p — sm a p 



sin a 



r 



cosa T 



(11) 



that makes the angle a p with the y axis directed along the polarization vector of the initial 
irradiation: E ini || y (see Fig. [p. 

Note that, in the limiting case of non- interacting molecules with U 2 = 0, the rotational 
mean-field FP equation §5§ reduces to the well-known linear rotational diffusion equation 
that has been widely used to study a large variety of problems based on the rotational 
diffusion model for the rotational motion of molecules in the presence of external fields. These 
include dielectric and Kerr effect relaxation of polar liquids (43-51|. rotational diffusion of 
a probe molecule dissolved in a liquid crystal phase [52-57] and the molecular reorientation 
in liquid crystal phases 58 6l| . 



2. Steady states 

The equilibrium angular distribution can generally be obtained as a stationary solution 
to the FP equation ()5]). It can be readily checked that the Boltzmann distribution 

/ st (u) = Z' 1 exp[-y(u)], Z st = J exp[-V(u)]du (12) 

determined by the effective potential gives the stationary solution. Note that the energy 
scale is determined by the temperature factor ksT incorporated into the rotation diffusion 
tensor, so that the internal energy, U[f), and the potential, V, are both dimensionless. 
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When the FP equation is linear, the stationary distribution ( TT2"j) describing the equi- 
librium state is unique. By contrast, this is no longer the case when the effective poten- 
tial ( imp depends on the elements of the averaged orientational order parameter tensor ([2]): 
V(u) = y(u|S). In this case, the components of the order parameter tensor in the stationary 
state, Sij = s!jj , can be found from the self-consistency equation 

Sij = J Qy(u)/ st (u|S)du. (13) 

which generally possesses several solutions representing multiple local extrema (stationary 
points) of the free energy 

F[f st ] = F st (S) = SySy - lnZ st (S). (14) 



Following the line of reasoning presented in Ref. 62] and using the effective free energy 
as the Lyapunov functional, it is not difficult to prove the H theorem for nonlinear FP 
equations of the form ([3]). It follows that, in the long time limit, all transient solutions 
converge to stationary ones. So, each stable stationary distribution is characterized by the 
basin of attraction giving orientational states (angular distributions) that evolve in time 
approaching the stationary distribution. 



B. Two-dimensional model of purely in-plane reorientation 

In the previous section our model has been formulated as the free energy FP equa- 
tion (jSJ) describing rotational diffusion of azo-dye molecules governed by the effective mean 
field potential (fTUj) . In this section, we concentrate on the limiting case of purely in-plane 
reorientation and deal with the simplified two-dimensional model 



Orf =d* <V + W = %f + o 3 %f V ) 



+ ffiV - Vd 2 J 



D (rat) t 



(15) 



can be derived by assuming that the out-of-plane component of the unit vector u de- 
scribing orientation of the azo-dye molecules is suppressed and, as is shown in Fig. [TJ, 
u = (sin0, cos0, 0). It implies that, similar to the model of a single axis rotator with 
two equivalent sites the molecules are constrained to be parallel to the substrate 

plane (the x-y plane). 



1. Potential 



When the out-of-plane reorientation is completely suppressed and u = (sin0, cos0, 0), 
the order parameter dyadic simplifies as follows 
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where Q2 is the two-dimensional in-plane part of the tensor. The in-plane order parameter 
tensor is given by the 2x2 block-matrix 

4Q 2 (0) = I 2 - 3 cos(20) a 3 + 3 sin(20) <n (17) 
expressed in terms of the Panli matrices: a 3 = (\ M and a x = (? j\ We can now 



v o -1 y 1 \ i io / 

substitute the expression for the order parameter tensor ([IE]) into the effective potential ( ITU]) 
and use the algebraic identity 

R r ot(-«)Q 2 (0)R r ot(a) = Q 2 (0 + a) (18) 

to derive the angular dependent part of the potential in the following form: 

V = (v* exp[2i<f)] + v exp[-2i<f>]) /2, (19a) 
v = Vi exp[— 2ia p ] + t> 2 (exp[2i0]), (19b) 
ui = 3w/Juv/4, v 2 = 9u 2 /8, (19c) 

/■27T 

where (...)= / ... /d0; v 1 and t> 2 are the photoexcitation and intermolecular interaction 
Jo 

parameters, respectively. 



Fourier harmonics and order parameters 

Our next step is to obtain the system of equations for the averaged complex harmonics, 
hn(T) — ( ex P[^0])( r )) that are proportional to the Fourier coefficients of the distribution 
function, /(</>, r) 

00 

/(0,r) = (27r)- 1 ^ /i„(r)exp[-m0]. (20) 

n=— 00 

To this end, we integrate the FP equation (|15|) multiplied by exp[m0] over the azimuthal 
angle. The resulting system reads 

d T h n (r) = -n 2 h n + n {v*(a p , h 2 )h n+2 {r) - v(a p , h 2 )h n ^ 2 (r)}, (21) 

where v(a p , h 2 ) = v = V\ exp[— 2ia p ] + v 2 h 2 . 

When /(</> + 7r) = /(</>), the odd numbered harmonics vanish, h 2k+1 = 0. For the even 
numbered harmonics, p k = h 2k , the system (12 ip can be conveniently recast into the form 

d rPk = -Ak 2 p k + 2k(v*p k+1 -i;p fc _i), A; = 1,2,... (22) 
Pk=Pk + m = (cos(2k(p)) + i{sm(2k(f))) = h 2k , p = 1, (23) 

where p\ = p\ + iq\ = (exp[2i0]) is the complex order parameter harmonics. 

For the in-plane order parameter tensor (fT7|) averaged over the azimuthal angle 

«W = S, = \ (* ;f Pl j + ?1 3pi ) = R«K) (" + °_) R„,(- Q ,), (24) 
*± = (l±3|pi|)/4, 2a e = aig(-pi*), (25) 
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the order parameter harmonics, pi, defines the principal values (eigenvalues), s + and s_, and 
the eigenvector, (cosa e , sina e ), for the largest eigenvalue s + giving the in-plane principal 
axis: n e = (cosa e , sina e , 0). So, the tensor (j24j) can be rewritten in the form of the in-plane 
order parameter dyadic: 

S 2 = s_I 2 + (s+ - s_)n e g> n e . (26) 

Temporal evolution of the in-plane order parameter tensor is determined by the order 
parameter harmonics, pi(t), that can be computed as a function of time by solving the 
initial value problem for the system (|22)) with the initial values of the harmonics Pfc(0) = pjj. 
representing the initial angular distribution /(</>, 0) = / m i (</>)• 



3. Symmetry 

It is not difficult to check that the system (|2ip is invariant under the symmetry transfor- 
mation 

h n — > h n = exp[ina]h n , v — >• v — exp[2ia]v = v(a p — a, h 2 ). (27) 

An important consequence of these relations concerns the special case where the reorienting 
light is switched off and V\ = 0. It can be readily seen that, at v — v 2 pi, a set of the 
harmonics, {pk}k, and that of the transformed harmonics {exp[2ika]pk}k both satisfy the 
system of equations (1221) . In particular, at Iuv — 0, the stationary angular distribution, 
f s t{4>), and the shifted (rotated) one, f s t{<f> — a), both represent steady states. 

By using the symmetry relations (1271) dynamical equations f[22|) can be conveniently 
changed into the system with the zero polarization azimuth, a v = 0. More precisely, from 
Eq. (127)1 . it follows that the transformed harmonics 

p k = exp[2ika p }p k (28) 

satisfy the system (1221) with the coupling coefficient v(a p ,p) and the initial distribution 
Zni(0) replaced by v = v(0, pi) and / n ; = / n i(0 — oc p ), respectively. Geometrically, such 
transformation can be described as a transition to the frame of reference where the y axis 
is directed along the polarization vector of reorienting light. 

So, in the subsequent section, we shall restrict our considerations to the special case of 
the system with the coupling coefficient 

v = vi + v 2 pi (29) 

taken at a p = 0. In this case, the polarization azimuth, a p , enters the initial data, p k = 
exp[2ika p ]p^ ni \ and the formula 

a e = a p + ^ arg(-p*) (30) 

gives the azimuthal angle of the principal axis ( 12 5 p expressed in terms of the order parameter 
harmonics p\. 
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Figure 2: (Color online) (a) Steady state order parameter harmonics p st as a function of 
the photoexcitation parameter v\ at t>2 = —2.5; (b) bifurcation curves in the v^-V\ plane 
are typical of the cusp catastrophe with the cusp singularity located at (—2, 0). 



III. KINETICS OF PHOTOINDUCED IN-PLANE REORDERING 



A. Bifurcations of equilibria 

The expression for the stationary distributions 

r-2-K 

f st = Z^expl-V] = Z^expHul cos(20 - a)], Z st = / exp[-V]d(j) (31) 

Jo 

comes from the general formula (I12p and represents the photosteady states in the two di- 
mensional case with the potential given in Eq. (j!9p where the coupling coefficient is changed 
from v to v — \v\ exp[za]. 

Equation ( 13 ID can now be combined with the relation 



exp[-wcos20] = I (v) + 2^2{-l) k I k {v)cos2k<f), (32) 

k=l 

where I k is the modified Bessel function of integer order [63], to derive the stationary state 
statistical integral 

Z st = 2ttI {\v\) (33) 
and the stationary state free energy (fH)) 



Fjv) 



\v\ 2 



+ Vi\v \ cos a 



-ln/ (H), (34) 
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where the additive constant is chosen so as to have the free energy vanishing at v = 0. 

The stationary points of the energy v = v st = v x + v 2 pf , define the steady state 
distributions (13"Tj) characterized by the stationary values of the order parameter harmonics, 
Pi = p{ +iqi \ which can also be found as solutions of the self-consistency equation ([TBI) . 
Close inspection of the formula flHj| shows that, for stationary points, the steady state 
coupling (potential strength) coefficient is real, {5 st =v± + v 2 p ( i t ' } =vi + v 2 p s t, with = 



and sina st = 0. So, we may closely follow the line of reasoning presented in Ref. [35j to 
perform the bifurcation analysis. 

According to Ref. [35], the steady state harmonics are given 

pi St) =pf = (-l) fc 4at)//ofe), (35) 
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where v st satisfies the self-consistency condition 

Pst = (vst ~ vi)/v2 = -Ii(v st )/I (v st ). (36) 

As is evident from the curve depicted in Fig. EJa), the number of solutions of the self- 
consistency equation (1361) varies between one and three depending on the values of the 
parameters v\ and v 2 . Figure EI a) shows the order parameter harmonics, = p st , plotted 
in the V\-p\ plane by using the parametrization 

and demonstrates that, for the intermolecular interaction parameter v 2 = —2.5 and suffi- 
ciently small photoexcitation parameters, the steady state free energy (|Hlj) has two local 
minima separated by the energy barrier and thus possesses three stationary points. More 
specifically, the steady states are determined by the free energy of the double-well potential 
form only if the inequalities 

v 2 < 4 c) = -2, «H < Vl < V W (38) 

are satisfied. 

The critical values of the parameter v\ depend on the reduced strength of intermolecular 
interaction v 2 and can be parametrized as follows 

vi =ui(f) = £-v 2 pi(Q, 

Geometrically, in the v 2 -Vi plane, equation (|39|) defines the bifurcation curves shown in 
Fig. [2](b). These curves form a bifurcation set which is the projection of the cusp surface 

= £-Cpi(0, 

Sb = { v 2 = C, (40) 

representing the bifurcation diagram in the three dimensional (v±,v 2 , pi) space (see Fig. [3D . 
Note that the cusp bifurcation occurs as a canonical model of a codimension 2 singularity [64( 



and the surface shown in Fig. [3] is typical of the cusp catastrophe [65j, [66 



B. Steady-state order parameter tensor 

From the discussion at the end of Sec. Ill B 31 and the relation (|35i) . at a v ^ 0, the 
stationary angular distributions are defined by the steady state harmonics 

P f ] = exp[-2tka p }(-l) k I k (d st )/I (d st ), (41) 

which are expressed in terms of the steady state coupling coefficient, v st — V\ + v 2 p s t 1 
that meets the self-consistency condition f[3"B"j) . For each steady state characterized by the 
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corresponding solution of Eq. (I3"6"j) . v st , the principal values of the in-plane order parameter 
tensor (J23J, s± t ' ) , can be computed from Eq. fTSol) by setting p 1 = p st = (v st — vi) / v 2 , whereas 
the azimuthal angle, ae , characterizing steady state orientation of the principal axis is 
given by the formula (13"0"j) with pi = p st : 

4 st) = (l±3|p st |)/4, (42a) 

(st) , 1 / \ J and ni st) 1 E uv , p st < 

«e '=« P + a axg(-p rt ) = <^ (st) . (42b) 

2 « p + 7T/2 and || E uv , p st > 




Irradiation time (min) 

(a) Order parameter harmonics (b) Angular distribution function 

Figure 4: (a) Order parameter harmonics, p\, computed as a function of dimensionless 
time parameter r by numerically solving the initial value problem for the system ( 12 2 p with 
the initial angular distribution ( l4"3"j) at a ini = 65 deg and p ini = —0.41. The photoexcitation 
and interaction parameters are: v i = 3 and t>2 = —1. The stationary state is unique and is 
characterized by the steady state harmonics p st — ~ —0.86. (b) Angular dependence 
of the orientational distribution function f l4"6]) . /(</>, r), at different values of the time 
parameter, r. The distribution is 7r-periodic and solid circles indicate the angles, 
0max = 7r /2 — « e , where peaks are located. 



An important conclusion to draw from this result is that, upon reaching the stationary 
state, the principal axis reorients along the direction perpendicular (parallel) to the po- 
larization vector of the reorienting light, e p = (— sincip, cosa p , 0), with a e — > ai st ' = a p 
(a e — > = 7r/2 + a p ) when the steady state order parameter harmonics is negative (pos- 
itive), p st = p[^ < (p st = p[^ > 0). Thus, in the photosaturated regime, the orientation 
of azo-dye molecules turned out to be defined solely by the polarization azimuth of the 
reorienting UV light. 

As it has been discussed in Sec. Ill B 3^ when the activating light is switched off and 
Iuv — (fi = 0), the steady states are continuously degenerate with the azimuthal angle 
of the principal axis playing the role of the degeneracy parameter. The reason for such 
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degeneracy is that the intermolecular interaction is taken in the rotationally invariant Maier- 
Saupe form The effect of the reorienting linearly polarized light is to lift the degeneracy 
by breaking the rotational symmetry of the intermolecular interaction. 




Figure 5: Angular dependence of the orientational distribution function, /, within the span 
of the period computed for various values of the time parameter, r = D rot t, at 
«ini = 65 deg and p; n i = —0.58. The photoexcitation and interaction parameters are 

v\ = 1.3 and t>2 = —5, respectively. 



C. Computational procedure and bistability effects 



According to the results of Ref. [38(, when the photoaligned azo-dye film is used as an 
aligning substrate in a NLC cell, orientation of the easy axis at the azo-dye substrate is 
dictated by the principal axis of the order parameter tensor f )24|) . In particular, under 
certain conditions, the difference between the angle a e and the easy axis azimuthal angle 
can be negligible. In this case, the kinetics of the easy axis can be modeled by solving the 
system f l2"2"j) and computing the angle a e (see Eq. flHUj) ) as a function of time. 

In our modeling, we numerically solve the initial value problem for the harmonics pk with 
the initial angular distribution 

/ ini oc exp[-w ini cos 2(0 - ini )]. (43) 

The distribution f )43p is characterized by the initial value of the order parameter harmonics 
Pi^ — Pmi = —h( v im) / Io( v mi) that gives the coupling coefficient v ini . At the initial instant 
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0.4 0.6 

Figure 6: Order parameter harmonics, pi, computed as a function of time parameter, 
t = D rot t, by numerically solving the initial value problem for the system 022]) with the 
initial angular distribution fl4"3"]) at a ini = 90 deg and a ini = 0.0 deg for — p ini = p u ± 0.05. 
The photoexcitation and interaction parameters are: v\ = 1.3 and V2 = —5. The stationary 
states are characterized by the steady state harmonics: p[ t ^ ~ —0.91, p^ ~ 0.77 and 

p u 0.48. 



of time, the harmonics are thus given by 

p fc (0) = exp[2i/ca ini ](-l) fc 4(t; ini )//oK 
In particular, for the order parameter harmonics, we have 



(44) 



Pi(0) = Pi n iCos(2a ini ), gi(0) = p ini sin(2a i: 



(45) 

The initial data defined in Eqs. f ]4"3"]) and (jUJ) represent the angular distribution of azo-dye 
molecules formed in the film after the preparatory stage of initial irradiation with linearly 
polarized UV light, E ini . Note that, when the polarization vector E ini is parallel to the y 
axis, the angle «i n i is equal to the polarization azimuth of the reorienting light a p . 

The system (122]) with the coupling coefficient v = v = v\ + v<if>\ describes how the 
harmonics evolve in time approaching, in the long time limit, the steady state values given 
in Eq. (135]) . These values are determined by the steady state order parameter harmonics 
that can be found by solving the self-consistency equation (136]) 

At sufficiently large positive photoexcitation parameter, V\ > \v c \, the equilibrium (photo- 
steady) state is unique and is characterized by the steady state order parameter harmonics 
p^ ■ This case is illustrated in figure Ha] that shows how the order parameter harmon- 
ics, pi = pi + iqi, relaxes to the equilibrium state, pi(t) — > p^ < and q\(t) — > at 
t — > oo. The equilibrium value of the principle (easy) axis angle ( 130]) then equals the 
polarization azimuth a p and a e (t) — > a p . The latter implies that, in the regime of pho- 
tosaturation, the azo-dye molecules align perpendicular to the polarization vector of the 
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reorienting, n e — > (cos a p , sin a p , 0) _L E^y. Such behavior is in complete agreement with 
our experimental data described in Sec. IIVI and has been previously observed in a number 
of experimental studies (see, e.g., the monograph 0] and references therein). 

Figure Ha] additionally illustrates that the magnitude of the order parameter harmonics, 
\pi\, related to the scalar order parameter, s + , given in Eq. (142 ap can be a nonmonotonic 
function of time. For the angular distribution function 



f(<P,r) = (2n)- 1 Pk(r)exp[-2zk(<p + a p )} 



(46) 



this implies that the distribution may experience broadening at the early stage of relaxation 
whereas it is getting narrower during the subsequent relaxation to the highly ordered photo- 
steady state. The effect of early stage disordering can be seen in Fig. [4b] that shows how the 
orientational distribution function (HBl) of azo-dye molecules evolves in time. Note, however, 
that such transient disordering is typically either less pronounced or completely suppressed. 
The latter case is demonstrated in Fig. [5] 

Temporal evolution of the angular distribution shown in Fig. is computed in the bista- 
bility region where the intermolecular interaction parameter, t> 2 , is below its critical value, 
t>2 = —2, and the photoexcitation parameter, V\ > 0, is sufficiently small, V\ < \v c \ (see 
Eq. f)38p ) . This is the region where, in addition to the equilibrium state with the negative 
order parameter harmonics p\ = p[ t \ the metastable state with the harmonics of opposite 
sign pi = p^ is formed. 



0.77 




Figure 7: Order parameter harmonics, pi, computed as a function of time parameter, 
t = D TOt t, by numerically solving the initial value problem for the system f ]22|) with the 
initial angular distribution ( l4U|) at = 90 deg and —pi n i = 0.65 for v\ = \v c \ ± 0.05. The 
interaction parameter and the critical value of the photoexcitation parameter are: v<i = —5 

and \v c \ ~ 1.51, respectively. 



Hence the stationary state free energy 
two local minima located at p^ and p^ 



takes the form of double-well potential with 
There is also a local maximum peaked at the 
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point p u in the interval ranged from to Referring to Fig. |5f a), at i>i > 0, this point 
with p u > lies on the upper half part of the branch representing unstable stationary states. 
In the case of real valued harmonics with = 0, it is the boundary point of the basins of 
attraction of the two steady states. In particular, this means that the order parameter 
harmonics p% converges to either p si ■* or p^ depending on whether its initial value pi(0) is 
below or above the boundary value p u . 

At a ini = and a ini = tt/2, the initial data (|44|) are represented by two sets of real valued 
harmonics. From Eq. (jUJ) it is evident that the only difference between the sets is the sign 
of the odd numbered harmonics . In particular, we have pi(0) = Pmi and pi(0) = — p ini at 
and a ini — tt/2, respectively. 



a: 



Otr. 



At a m i = and p- m i < 0, the harmonics characterizing the equilibrium state, p[ t \ is 
the attracting point for the order parameter harmonics and the azimuthal angle of the 
principal (easy) axis, remains intact, ai st ' = 0. The latter comes as no surprise because 
the polarization vectors of the initial irradiation and of the reorienting light are parallel, 

EjjV || Ejni 

Interestingly, orientation of the easy axis may stay unchanged even if the vectors, E^y 

-n/2. This occurs when the initial ordering is 
> p u . As demonstrated in Fig. [6l in this case, 



and E ini , are perpendicular and = a p 
high, so that — p irii is greater than p u , —p h 

the attracting point is the metastable harmonics p^ rather than the equilibrium one 
and «l st ' = 0. 
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Figure 8: Order parameter harmonics, pi, computed as a function of time parameter, 
t = D rot t, at a ini = 89.8 deg and p ini = — pj£' . Other parameters are listed in the caption 

of Fig. El 



It can also be seen from the curves of Fig. [6] that, in the opposite case of low ordered initial 
states with —p m i < p u , the order parameter harmonics p\ relaxes back to the equilibrium 
value p^' '. So, at —p m \ < Pu-, the azimuthal angle, = a p = tt/2, defines the steady state 
easy axis which is normal to the direction of its initial orientation. 
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In Fig. [7J we show that, when the initial ordering is sufficiently high, the process of 
reorientation crucially depends on whether the photoexcitation parameter, V\, which is pro- 
portional to the intensity of the reorienting light, Iuv, exceeds its critical value, \v c \. This 
critical value determines the intensity threshold, I c = 4|t> c /(3u/)|, and the easy axis reori- 
entation is suppressed at a p = it/ 2 in the low intensity region below the threshold where 
Iuv < Ic- 

Similar threshold behavior has been previously observed in a number of experimental 
studies on light-induced easy axis reorientation in nematic liquid crystal cells with photo- 
sensitive poly mer substrates [3, 0, 68|. For azo-dye films, the experimental data of 



Refs. [42, |69j also suggest the presence of a threshold when the polarization vector E^y is 
directed along the initial easy axis (see the subsequent section for more details). Our result 
that under certain conditions highly ordered initial distributions may relax to the metastable 
stationary state offers a possible explanation for this. 

Referring to Fig. [HI it is seen that small deviations of the angle a- m \ from it/2 have a 
destructive effect on metastability of the attracting point even if — p ini = p^ . It turned out 
that the basin of attraction of the metastable state is unstable with respect to fluctuations 
of the polarization azimuth a p . In the vicinity of the critical angle a- m i = tt/2, the presence 
of the metastable state manifests itself as appreciably retarded relaxation to the equilibrium 
state (see Fig. [HJ. Similarly, the shape of the azimuthal angle (130|) vs time curve shown in 
Fig. M points to a time delay in the process of principal (easy) axis reorientation. 




Figure 9: Principal axis azimuthal angle, a e , computed as a function of time parameter, 
r = D Iot t, at a ini = a p = 89.8 deg and p ini = — . Other parameters are listed in the 

caption of Fig. 
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IV. THEORY VERSUS EXPERIMENT 



In the previous section we have examined the experimentally important predictions of 
the model. These predictions can, in principle, be tested and, in this section, we discuss 
some of our recent experimental results on the reorient at ional dynamics of the electrically 
assisted light-induced azimuthal gliding of the easy axis at varying polarization azimuth of 
LPUV reorienting light j42|. 

Such gliding is found to take place on photoaligned azo-dye layers when irradiation of 
nematic LC (NLC) cells with LPUV light is combined with the application of ac in-plane 



electric field |69|, |70j. Since the effects of electric field is well beyond the scope of the 
theoretical approach under consideration, the regime of purely light-induced reorientation 
will be our primary concern and we deal with the data measured in the electric field-free 
regions 

In our experiments, LC cells (d = 17 A ± 0.2 jttm) of sandwich like type were assembled 
between two amorphous glass plates. The upper glass plate was covered with a rubbed 
polyimide film giving the strong planar anchoring conditions. A film of the azobenzene 
sulfric dye SD1 (Dainippon Ink and Chemicals) js[ was deposited onto the bottom substrate 
on which transparent indium tin oxide (ITO) electrodes were placed. The electrodes and 
the interelectrode stripes (the gap was about g = 50 /iin and the in-plane ac voltage was 
U = 100 V with / = 3 kHz) were arranged to be parallel to the direction of rubbing (the x 
axis). 

The azo-dye SD1 layer was initially illuminated by LPUV light at the wavelength A = 
365 nm. The preliminary irradiation produced the zones of different energy dose exposure 
D p = 0.27, 0.55 J/cm 2 characterized by relatively weak azimuthal anchoring strength. The 
light propagating along the normal to the substrates (the z axis) was selected by an interfer- 
ence filter. Initial orientation of the polarization vector of the actinic light, E ini , was chosen 
so as to align azo-dye molecules at a small angle of 4 degrees to the x axis, ipo 4 deg. 

The LC cell was filled with the nematic LC mixture E7 (Merck) in isotropic phase and 
then slowly cooled down to room temperature. Thus we prepared the LC cell with a weakly 
twisted planar orientational structure where the director at the bottom surface no is clock- 
wise rotated through the initial twist angle tpo ~ 4 deg which is the angle between no and 
the director at the upper substrate (the 

In addition to the electric field, E = 2 V/fim, the cell was irradiated with the reorienting 
LPUV light beam (Iuv — 0.26 mW/cm 2 and A = 365 nm) normally impinging onto the 
bottom substrate. For this secondary LPUV irradiation, orientation of the polarization 
plane is determined by the polarization azimuth, a p . 



Our experimental method has already been described in Refs. [69l . 170]. In this method, 
NLC orientational structures were observed via a polarizing microscope connected with a 
digital camera and a fiber optics spectrometer. The rotating polarizer technique was used 
to measure the azimuthal angle (p e characterizing orientation of the easy axis. In order to 
register microscopic images and to measure the value of (p e , the electric field and the reori- 
enting light were switched off for about 1 min. This time interval is short enough to ensure 
that orientation of the easy axis remains essentially intact in the course of measurements. 
The measurements were carried out at a temperature of 26°C. 

The case where the reorienting light is linearly polarized along the initial surface director 
n and a p = tt/2, was studied in Refs. (69|, 70[. It turned out that, by contrast to the 



electrically assisted light-induced gliding, the purely photoinduced reorientation is almost 
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Irradiation time (min) 

Figure 10: Easy axis angle versus the irradiation time at different values of the 
polarization azimuthal angle of the reorienting UV light a p . The experimental data were 
measured by using the linearly polarized activating light with A = 365 nm and 
/ = 2.6 W/m 2 . Solid lines represent the theoretical curves computed by numerically 
solving the initial value problem for the system f[2"2"l) with the initial angular 
distribution (T4"3l . The parameters used in calculations are as follows: the photoexcitation 
(interaction) parameter is v\ = 1.3 (v 2 = —5); p- mi = p[^ (a p = 25 deg and a p = 43 deg) 
and p ini = —0.25 at a p = 65 deg. Other parameters are listed in the caption of Fig. EJ 



entirely inhibited at a p = it/2. As it can be seen from Fig. [TU], the latter is no longer the 
case for the reorienting light with a p ^ ir/2. 

Thus the result is that, as opposed to the cases with non- vanishing electric field and 
a p 7^ 7r/2, the purely photoinduced reorientation of the easy axis has been completely 
suppressed when the polarization plane is parallel to the initial easy axis. From the above 
discussion, our model predicts similar behavior that occurs in the bistability region where 
the reorientation is characterized by the threshold for the intensity of reorienting light. 

Figure [10] shows the easy axis angle measured as a function of the irradiation time at 
various values of the polarization azimuth. The measurements were performed in the region 
outside the interelectrode gaps where the ac electric field is negligibly small. In the zero- 
field curves, the easy axis angle increases with the irradiation time starting from the angle of 
initial twist, (po, and approaches the photosteady state characterized by the photosaturated 
value of the angle close to a p . 

The experimental data can be fitted by the theoretical curves computed from the for- 
mula (1301) . where the order parameter harmonics, p\, is found by numerically solving the 
initial value problem for the system fl22l) with the initial angular distribution (j43l . In Fig. [10] 
we make a comparison between the experimental data and the theoretical curves computed 
in the bistability region by setting the photoexcitation and interaction parameters, v\ and 
i>2) equal to 1.3 and —5, respectively. In this region, the equilibrium state is character- 
ized by the steady state harmonics p st = ~ —0.91, whereas the stationary harmonics 
m 0.77 corresponds to the metastable state. For a p = 25 deg and a p = 43 deg, the 
curves were computed at p ini = p^ . It follows that |pi(0)| = |p st | and the angle a ini is as- 
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sumed to be the only parameter describing difference between the initial and the equilibrium 
angular distributions. The rotational diffusion coefficient D TOt can be estimated at about 
1.6 x 10~ 4 s _1 Using this value to fit the data measured at a p = 65 deg, we find that the 
theoretical and experimental results are in agreement when the initial ordering parameter 
Pini ~ —0.25 differs from the equilibrium value. 
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Figure 11: Easy axis angle versus the irradiation time in the electric field- free region at 
different values of the polarization azimuthal angle of the reorienting UV light. The 
experimental data were measured by using the linearly polarized activating light with 
A = 365 nm and / = 2.6 W/m 2 . Solid lines represent the theoretical curves computed from 
the formula (1481 at t e = 7 e /W / e = 23 min (for ip m = 25 deg and ip m = 43 deg), and 

t e pa 14.4 min. for ip m = 65 deg. 
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We conclude this section with the remark concerning the phenomenological model formu- 



lated in Refs. [41|, |70j to describe the effect of the electrically assisted photoinduced gliding of 
the easy axis. According to this model, in the absence of ac electric field, temporal evolution 
of the easy axis azimuthal angle ip e can be described by the balance torque equation [42] 

le-^- = — — sm2(^ e -<p m ), (47) 

where 7 e is the specific viscosity of gliding; W e is the effective anchoring parameter which 
defines the strength of coupling between the easy axis n e and the initial state of surface 
orientation; ip m is the polarization dependent phase shift. The solution of the dynamic 
equation (1471) 

tan(^ e (t) - ip m ) = tan(^ e (0) - ip m ) exp[-t/t e ], (48) 

where t e = , y e /W e the characteristic time of easy axis reorientation, can now be used to fit 
the experimental data. The results are presented in Fig. ITTl 

From the curves depicted in Fig. [11] it can be inferred that the phase shift ip m is close 
to the polarization azimuth of the reorienting light. The value of the characteristic time, 
t e pa 14.4 min, estimated for (p m = 65 deg differs from the one t e ~ 23 min obtained for other 
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two angles, ip m = 25 deg and <p m = 43 deg. From the other hand, in Fig. [TUJ the case with 
a p ~ 65 deg is distinctive in initial ordering of the film. This result points to the fact that the 
coupling strength of the easy axis W e depends on the photoinduced order parameter of the 
initially irradiated azo-dye layer. Such dependence was one of the key assumptions taken in 
Ref. [&9[ when applying the phenomenological model jU 70[ to interpret the experimental 
data for photoaligned azo-dye layers prepared at different initial irradiation doses. 



V. DISCUSSION AND CONCLUSIONS 

In this paper we have studied the kinetics of orientational reordering that takes place in 
the azo-dye films under the action of LPUV reorienting light. In our theoretical approach, 
the process of reordering is described in terms of a rotational Brownian motion and we have 
employed the 2D diffusion model formulated in Ref. [35] to explore the peculiarities of the 
photoinduced reordering in the regime of purely in-plane reorientation. 

Orientational ordering of azo-dye molecules is characterized by the in-plane (2D) order 
parameter tensor f[24"j) . This tensor is expressed in terms of the order parameter harmonics, 
pi, that enters the formula fl25|) giving both the value of the scalar order parameter, s+, and 
the azimuthal angle, a e , of the principal (easy) axis, n e . 

It is found that dynamics of the order parameter harmonics essentially depends on the 
initial ordering and the polarization azimuth of the LP reorienting light. In the long time 
limit, the order parameter harmonics approaches its steady state value reaching the regime of 
photosaturation. We have proved that, in the steady state, the easy axis azimuthal angle is 
dictated solely by the polarization azimuth, Eq. (1421) ) . whereas the photostationary 

value of the scalar order parameter is determined by the photoexcitation and intermolecular 
interaction parameters, V\ oc and t> 2 , defined in Eq. fl L9c|) . 

In the t>i-t>2 plane, inequalities ( 13"5|) define the bistability region (see Fig. I^b)) with 
two locally stable stationary states: equilibrium and metastable ones. An experimentally 
important bistability effect is that the easy axis reorientation is inhibited when the intensity 
of the reorienting light linearly polarization along the initial easy axis is below its critical 
(threshold) value. Such threshold behavior was observed in our experiments where the easy 
axis angle was measured as a function of the irradiation time at varying polarization azimuth 
of LPUV reorienting light. 

An alternative phenomenological model that also predicts the presence of the reorienta- 



tion threshold was suggested in Ref. [21|. In that model, the threshold effect comes about 
from the interplay between two competing easy axes. By contrast, in our model, the bistabil- 
ity appears as an intrinsic property of the reordering process in azo-dye films that manifests 
itself as the threshold of reorientation. 

There is another interesting bistability effect that occurs when the light is switched off 
(Iuv — 0) and t>2 < = — 2. In this case, owing to the rotational symmetry of the two- 
particle interaction ([9]), the easy axis angle, a e , becomes undetermined and plays the role of 
the degeneracy parameter. The LPUV light field lifts the degeneracy and the dynamics of 
the angle turned out to be predominately governed by the photoexcitation parameter, v\. 

For the scalar order parameter which is related to the magnitude of the order parameter 
harmonics, \pi\, the rate of relaxation is typically at least twice faster than that for the 
easy axis angle. As is illustrated in Fig. HaJ depending on the initial conditions, temporal 
evolution of this order parameter can be non-monotonic. These effects can also be seen from 
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Figs. |4b] and that show how the orientational angular distribution of azo-dye molecules 
evolves in time. 

By numerically solving the dynamical system for the angular distribution harmonics (122]) 
we have fitted the experimental data on the light-induced azimuthal gliding of the easy 
axis at the photoaligned azo-dye SD1 layer. In our analysis, the photoinduced reordering of 
azo-dye molecules was assumed to be the key factor dominating the easy axis reorientation. 

Despite the fact that good agreement between the computed curves and the data (see 
Fig. IT0|) counts in favor of our simplified approach, a deeper insight into complicated molec- 
ular mechanisms behind the gliding under consideration requires additional experimental 
studies and a more systematic theoretical treatment of the processes involved. Such treat- 
ment has to deal with the interplay of photoinduced ordering in azo-dye films, the anchoring 
energy effects [38[ and the adsorption-desorption processes underlying the gliding 

phenomenon. 
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Potential 2.1.1/5873; Grant NK-410P; HKUST CERG RPC07/08.EG01 and CERG 612208, 
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